Complete dynamics of H2 in strong laser fields 
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Based on a combined quantum-classical treatment, a complete study of the strong field dynamics 
of H^", i.e. including all nuclear and electronic DOF as well as dissociation and ionization, is 
presented. We find that the ro- vibrational nuclear dynamics enhances dissociation and, at the same 
time, suppresses ionization, confirming experimental observations by I. Ben-Itzhak et al. [Phys. 
Rev. Lett. 95, 073002 (2005)]. In addition and counter-intuitively, it is shown that for large initial 
vibrational excitation ionization takes place favorably at large angles between the laser polarization 
and molecular axis. A local ionization model delivers a transparent explanation of these findings. 



Even for nature's simplest molecule H^", no full- 
dimensional description (including rotation of the nu- 
clei) of the ionization dynamics has been given up to 
date. While in the early days important mechanisms 
of strong field ionization could be identified in calcula- 
tions with frozen nuclei [![, it was shown subsequently 
that nuclear vibration does have a crucial influence on 
the ionization dynamics. Prominent effects include the 
washout of charge resonance enhanced ionization at spe- 
cific internuclear distances (CREI) 0, Q by vibrational 
motion 0, [f| as well as Lochfrafi (internuclear separa- 
tion dependent ionization) 0, Q • For laser pulses much 
shorter than a typical rotational period it has been rea- 
sonably assumed that rotation can be neglected 0, Q. 
However, recent theoretical work has revealed the 
importance of rotation for the dissociation channel even 
for very short laser pulses. From a key experiment [11 1 
it could be concluded that rotation and/or orientation is 
also relevant for ionization, which appears to be overesti- 
mated when the molecule remains aligned with the laser 
polarization. 

All these findings can be put into a stringent per- 
spective with the results from our present calculations 
taking all degrees of freedom into account. The mixed 
quantum-classical non-adiabatic quantum molecular dy- 
namics (NA-QMD) method (i1j-fl4| which has enabled 
us to perform a full dimensional calculation is still ap- 
proximate since it treats the nuclei classically. Hence, we 
have carefully assessed its validity for parameter regimes 
where accurate quantum calculations can be performed. 
In [9| dissociation and ionization of laser-aligned Hj" 
molecules were calculated and in we have obtained 
vibrationally resolved angular distributions of dissociated 
fragments for laser parameters where ionization can be 
neglected. In both cases we have found quantitative 
agreement with appropriate quantum results, in the first 
case from a dimensionally reduced quantum approach Q , 
in the second case from full-dimensional quantum calcu- 
lations, however restricted to dissociation [l4j. There- 
fore, the mixed quantum-classical NA-QMD approach 
appears to be well suited to quantify ionization for H^" 



including all degrees of freedom, for which no calculations 
exist so far. 

Our starting point is the time-dependent electronic 
Schrodingcr equation (atomic units are used) 



i^*(r,*;{R-j}) 



A=l 



<f>(r,t;{R A }), (f) 



which is solved in basis expansion sclf-consistcntly with 
Newton's equations of motion for the nuclei. Here, 
(A = 1, 2) and r denote the nuclear and electronic coor- 
dinates, T T = — 1/2V 2 the electronic kinetic energy, and 
e(t) is the electric field of the laser. Details regarding the 
method and its numerical implementation can be found 
in [3, EH ■ At the final time tf the vibrationally resolved 
ionization (I) and dissociation (D) probabilities are con- 
structed from probability densities along relevant nuclear 
trajectories, 



P l (y) = l-Z- 1 Y^N\R i ](t t ) 



(2) 



1=1 



P D (u)=Z- 1 J2 ^K^f), (3) 

i:Ri(t { )>Rr, 

where AT[R](t) is the norm of the electronic wavefunc- 
tion, which is a functional of the trajectory R(t) and a 
function of time. Its decrease in time, due to the pres- 
ence of an absorber potential [lU 14 1, is a measure of 



ionization. Trajectories with final internuclear distance 
R > i?D = f a.u. contribute to dissociation. Z v denotes 
the number of trajectories, sampled according to a micro- 
canonical distribution of a given initial vibrational state 
v. With the molecule initially in its rotational ground 
state, the nuclear rotation angle 9 is sampled uniformly 
in the interval [0, where one has to take into account 
the additional weighting factor sin 6i (to) [lij ]. Comple- 
mentary to the full dimensional calculations we have also 
performed dynamical calculations but with fixed nuclear 
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orientation (frozen rotation) in order to study the influ- 
ence of the rotational dynamics on the results. 

We have determined vibrationally resolved ionization 
and dissociation probabilities for a linearly polarized laser 
pulse with an amplitude e(i) = eo sin 2 (nt/T) cos uit at a 
wavelength of 800 nm (corresponding to u w 0.057 a.u.), 
peak intensity of J = 2 • 10 14 ^r (corresponding to eo w 
0.075 a.u.) and total pulse length of T = 50 fs. Similar 



laser parameters were used in recent experiments llj, |15 , 

Fig.QJi reveals that the full dimensional integral results 
for ionization (red circles, including rotational dynam- 
ics) agree quite well with those for frozen nuclear axis 
(dashed lines), where in both calculations the molecules 
are isotropically distributed in the beginning. This is 
understandable since over the length of the laser pulse 
(50 fs) during which ionization is possible the molecule 
hardly rotates. Dissociation, however, is much slower and 
its probability consequently deviates considerably from 
the result for frozen angles (black squares versus black 
dashed line). Only for weak binding at large v dissocia- 
tion happens even more quickly than ionization so that 
Pd agrees with the frozen result. Moreover, this implies 
that beyond fragmentation saturation (Pi + Pd = 1, ver- 
tical dotted line) , ionization is suppressed by dissociation 
according to Pi(Pd) = 1— Pd- Finally, the two maxima in 
Pd around v = 4 and v = 11 above the frozen probability 
are due to the well known dressed two- and one-photon 
states, respectively |14|. 

The corresponding probabilities for the molecule 
aligned with the laser polarization (Fig. [T]d) differ quan- 
titatively and qualitatively from the full dimensional 
analogs (Fig. |TJl) . Due to the much stronger ioniza- 
tion, fragmentation saturation sets in already near v = 4. 
Hence, for the aligned molecule ionization suppresses dis- 
sociation for high v : according to Pd(Pi) = 1 — Pi, in 
striking contrast to the real case (Fig. QJi). For the high- 
est vibrational levels ionization decreases again, because 
there is a large probability that the molecules have an 
internuclear distance outside the strong ionization region 
(Fig. [2]) when the laser reaches peak intensity. From 
Fig. [2] it is also clear why the full dimensional and laser 
aligned molecular response to the light pulse is so dif- 
ferent: The coupling to the light changes considerably 
for different alignment angles and it is strongest in the 
aligned case where molecular enhancement mechanisms 
such as dressed state resonances and CREI are operative. 

The angular distributions of dissociated and ionized 
fragments (Fig. [3} confirm what we have concluded from 
the integral probabilities of Fig. [1] For the ionization 
channel, the angular distributions follow roughly the be- 
havior expected from geometric alignment. In contrast, 
the angular distributions for dissociation from full dimen- 
sional and rotationally frozen calculations exhibit an op- 
posite trend for increasing angle. Moreover, the strong 
alignment (for low and high v) and anti-alignment for in- 
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FIG. 1: (Color online) Ionization probability Pi(v) (red cir- 
cles), and dissociation probability Po(y) (black squares) for 
(a) initially rotationally cold molecules (isotropically orien- 
tated) and (b) laser-aligned molecules (see also Q|) as a func- 
tion of the initial vibrational state v. The dashed lines in 
(a) represent results for fixed nuclear orientation (frozen ro- 
tation). The vertical line denotes the vibrational level, above 
which fragmentation saturates. Pd + Pi = 1. Finally, the ion- 
ization probability based on the LIM (Eq. [6]) is shown with 
blue triangles. 
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FIG. 2: (Color online) The ionization probability P(R) = 
1 — iV(R, T) calculated for fixed nuclear positions with the 
parallel and perpendicular component Ru — Rcos8 resp. 
R± — Psin#. Geometric alignment roughly manifests itself in 
increasing ionization probabilities towards the laser polariza- 
tion direction. Near laser polarization direction clearly a re- 
gion of strong ionization is visible, in particular the 1-photon- 
resonance at R ~ 4.3 a.u. as well as additional enhanced ion- 
ization maxima related to CREI. 



termediatc v which has been explained and understood 
within the Floquct picture (e.g., |l4j) is missing for frozen 
angles. 

Summarizing the interpretation of the results so far 
it has become clear that a full dimensional approach is 
necessary to capture the dynamics adequately. Neverthe- 
less, as the ionization probability depends exponentially 
on the laser intensity, a simplified dynamical description 
may be possible. In our case Pi from Fig. [2] is reduced 
ten times when using half the peak intensity which is also 
observed experimentally [llj]. Hence, the full dynamical 
ionization should mainly be determined by the nuclear 
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FIG. 3: (Color online) Angular distribution of fragments 
starting from vibrational levels v = 4, 9, and 14 for ion- 
ization -Pi(#) (red bars) and dissociation Pd(0) (green bars) 
from full dimensional calculations (cf. [TJ] for definition of 
Pd(0), Pi(0) analogous). The results calculated with frozen 
rotation are given with gray circles. 



positions around peak intensity (25 fs into the pulse), 
dependent on the initial vibrational state as illustrated 
in Fig. 0^. To understand the influence of nuclear motion 
on ionization, we define the ionized nuclear density 



pi(R-) 



i=i n 



dt JV[Ri](t) <5(R-R(i)) , (4) 



which reveals preferential nuclear positions for ionization 
(Fig. |3)d). Surprisingly, for large v ionization proceeds 
favorably at large angles. It also corroborates the "naive" 
picture of Fig. |4^i that ionization mainly takes place at 
the time of peak intensity and is therefore sensitive to 
the dynamically reached nuclear positions at this time. 

To quantify this observation, we may assume that ion- 
ization happens instantaneously neglecting memory ef- 
fects for nuclear motion in the approximation 



iV[R i ](i)«-r(R J (t),t) N[Ri](t) 



(5) 



with an instantaneous local ionization "rate" r(R, t) = 
— iV(R, t) /N(R, t). Note that this assumption is not triv- 
ial, since the ionization of any trajectory in the swarm 
depends in principle on its whole time evolution. The 
instantaneous local ionization rate r(R, t) is calculated 
within the fixed-nuclei approximation using the same 
laser conditions as in the full dynamical calculations. 
This rate is not a simple function of time, internuclear 
distance and alignment angle, but also depends on the 
pulse shape of the laser. Inserting §5§ into (@| leads to a 
local ionization model (LIM), where the ionized nuclear 




FIG. 4: (Color online) (a) Trajectories of fragmenting 
molecules (red lines) as well as nuclear positions at peak 
intensity (black dots) for exemplary initial vibrational lev- 
els (y — 4,9,14). (b) Corresponding exact ionized density 
pi(R) and (c) ionized density in the LIM pp IM (R) (linear color 
scale). The nuclear coordinates are the same as in Fig. [2] 



density takes the product form 



p x LIM (R) = J dtr(Rt) p(R,t) 



(0) 



with the time-dependent nuclear density 



p(R,t) = Z-^NiRijit) 5(R — Ri(t)) . (7) 

i=l 

As one can see from Fig. \$jp, the separation of nuclear 
dynamics and ionization according to (|6]) reproduces the 
full result Fig. quite well and is in accord with the 
spirit of [rjj with the difference that our ionization rates 
are extracted from the exact laser pulse. 

Note, however, that the LIM overestimates ionization 
for low v, which is due to the proximity of a 1-photon res- 
onance close to R « 4.3 a.u. For higher v the model re- 
produces that ionization takes place dominantly at inter- 
mediate angles as a consequence of dissociation dynam- 
ics [14|. The quantitative accuracy of the LIM (Fig. QJi, 



blue triangles) is overall comparable to the rotationally 
frozen result (dashed line), however LIM better reflects 
dynamical details, such as the influence of dressed state 
resonances. 

A serious quantitative comparison with the experi- 
ment ll|, which reveals one order of magnitude differ- 
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FIG. 5: (Color online) Franck- Condon averaged angular dis- 
tribution for ionization (Pi(9)) and dissociation (Pd(6)). 



ence between total ionization and dissociation yield, re- 
quires Franck Condon averaging over initial vibrational 
levels and in addition focal volume and thermal averag- 
ing. While the former is a generally valid result, the lat- 
ter depends on the exact conditions of individual experi- 
ments, which differ and are not trivial to quantify. Hence, 
we restrict ourselves to calculate full Franck-Condon av- 
eraged [l~8| angular distributions for ionization and dis- 
sociation and discuss modifications due to volume aver- 
aging qualitatively. As one can see in Fig. [5j Pd(#) is 
more strongly aligned along the laser polarization than 
P\{9) in seeming discrepancy with the opposite obser- 
vation in the experiment. However, as noted in [llj ]. 
the dynamics changes substantially with small changes 
in pulse length if the latter is close to the vibrational 
time scale which is the case here. Hence, the small dif- 
ferences in experimental and theoretical pulse length and 
shape may matter. Secondly, the spatial intensity profile 
of the laser will favor contributions of ionized fragments 
in the wings of the laser focus which originate dominantly 
from laser-aligned molecules. In the experiment [U| low 
laser intensities contribute considerably to the signal due 
to the width of the ion beam used. On the other hand, 
Pd(#) would broaden for large contributions from lower 
intensities (cf. the IDS procedure used in [3]). The 
total dissociation and ionization probabilities from our 
calculation, Pi = 0.13 and Pd = 0.41, reproduce the ex- 
perimental trend, that ionization is considerably smaller 
than dissociation in sharp contrast to predictions from 
dimensionally reduced calculations. 

In summary, we have presented a complete study of 
strong field ionization and dissociation of H \ . All nuclear 
and electronic degrees of freedom have been included in 
the framework of the mixed quantum-classical NA-QMD 
method whose applicability has been checked carefully 
in prior work by comparison to quantum results where 
available. It is worthwhile to note, that NA-QMD fl2l.HH 
is also applicable to polyatomic many-electron systems 

We have found, that nuclear rotation enhances disso- 
ciation so that Pd > Pi for H^ under short pulses in 



agreement with the experiment [ill ] but in discrepancy 
with dimensionally reduced calculations. Furthermore, 
vibrationally resolved ionization P\(y) is reasonably well 
reproduced with frozen rotation, while Pd (y) is enhanced 
in comparison to frozen nuclear geometry with two max- 
ima whose origin can be easily traced to one- and two- 
photon dressed states. The latter connection is almost 
completely lost in the dimensionally reduced result since 
due to the strong ionization, fragmentation saturation 
occurs already for v > 4 and leads to suppression of dis- 
sociation. 

We conclude that full dimensional calculations are nec- 
essary for short pulses to obtain even qualitatively correct 
results. Nevertheless, a simplification arises in a local 
ionization model (LIM) from the separation of ioniza- 
tion and nuclear dynamics, taking advantage of the fact 
that ionization is fast, but happens dominantly at peak 
intensity, when nuclei have already moved. Finally, we 
could demonstrate that rotationally frozen dynamics is a 
good approximation when dissociation or ionization hap- 
pens fast with respect to the rotational time scale. This 
is always the case for ionization with a 50 fs pulse and 
applies regarding dissociation to the highest vibrational 
levels, which dissociate before the molecule rotates sub- 
stantially. 
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